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PREFACE 


This  investigation  was  conducted  at  the  U.  S.  Army  Engineer  Water- 
ways Experiment  Station  (WES)  under  the  sponsorship  of  the  Office,  Chief 
of  Engineers,  Department  of  the  Army,  as  a part  of  Project  CWIS  311^5, 
"Liquefaction  Potential  of  Dams  and  Foundations  Under  Earthquake 
Excitation." 

The  investigation  was  conducted  by  Drs.  G.  Y.  Baladi  and  B.  Rohani 
during  the  period  September  1976  through  August  1977,  under  the  general 
direction  of  Mr.  J.  P.  Sale,  Chief,  Soils  and  Pavements  Laboratory,  and 
Drs.  J.  G.  Jackson,  Jr.,  Chief,  Soil  Dynamics  Division,  F.  G.  McLean, 
Chief,  Earthquake  Engineering  and  Vibrations  Division  (EE&VD),  and 
W.  F.  Marcuson  III,  Project  Leader  (EE&VD).  The  report  was  written 
by  Drs.  Baladi  and  Rohani. 

COL  J.  L.  Cannon,  CE,  was  Director  of  WES  during  the  preparation 
of  this  report.  Technical  Director  was  Mr-;'V-.  R.  Brown. 
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LIQUEFACTION  POTENTIAL  OF  DAMS  AND  FOUNDATIONS 


DEVELOPMENT  OF  A CONSTITUTIVE  RELATION  FOR  SIMULATING 
THE  RESPONSE  OF  SATURATED  COHESIONLESS  SOIL 

PART  I : INTRODUCTION 

Backgroiuid 

1.  During  1975  a research  effort  was  initiated  at  the  U.  S. 

Army  Engineer  Waterways  Experiment  Station  (WES)  to  develop  an  analyti- 
cal tool  for  studying  the  liquefaction  potential  of  dams  and  foundations 
subjected  to  earthquake  excitation.  Specifically,  the  objective  was  to 
develop  a rational  computational  method  (based  on  continuum  mechanics 
precepts)  for  the  performance  of  effective  stress  analyses  of  realisti- 
cally posed  bovmdary-value  problems  involving  fxilly  saturated  cohesion- 
less soils.* 

2.  The  investigation  was  divided  into  three  phases.  Phase  I, 
completed  in  1976,^  involved  the  development  of  a three-dimensional 
(3D)  elastic-plastic  isotropic  constitutive  model  that  qualitatively 
simulated  some  of  the  basic  stress-strain-pore  pressure  response  fea- 
tures observed  in  laboratory  tests  conducted  on  fully  saturated  cohe- 
sionless soils.  This  model,  however,  does  not  treat  observed  strain- 
softening behavior,  nor  does  it  predict  the  progressive  increases  in 
pore  pressure  observed  under  low-amplitude  (subfailure)  cyclic  shear 
loading  conditions — two  phenomena  considered  important  contributors  to 
the  liquefaction  problems.  These  two  features  were  successfully  modeled 
during  Phase  II  of  the  investigation  by  appropriate  extensions  of  the 
Phase  I constitutive  model.  The  objective  of  Phase  III  (not  yet  begun) 
is  to  incorporate  the  Phase  II  version  of  the  constitutive  model  into  a 
suitable  numerical  code  for  future  use  in  the  conducting  of 

* The  terms  "cohesionless  soil"  and  "sand"  are  used  interchangeably 
throughout  this  report  to  designate  earth  materials  with  negligible 
cohesion. 
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effective-stress  analyses  and  assessment  of  liquefaction  potential  of 
various  earth  structures  subjected  to  earthquake  or  other  dynamic 
loading  conditions. 


Purpose 

3.  The  purpose  of  this  report  is  to  document  the  results  of 
Phase  II  of  the  investigation,  i.e.,  the  extension  of  the  Phase  I model 
to  (a)  include  strain-softening  shear  behavior  and  (b)  predict  the 
progressive  increase  of  pore  pressure  observed  when  a saturated  sand  is 
subjected  to  undrained  cyclic  shear  loadings. 

Scope 

k.  The  typical  mechanical  behavior  of  saturated  coehsionless 
soils  subjected  to  various  laboratory  test  loading  conditions  is  dis- 
cussed in  Part  II.  Part  III  presents  the  mathematical  formulation  of 
the  constitutive  model.  The  behavior  of  the  constitutive  model  under 
simulated  triaxial  test  conditions  is  demonstrated  in  Part  IV.  Part  V 
presents  comparisons  of  model  predictions  with  the  measured  laboratory 
response  of  two  sands.  A summary  and  recommendations  are  given  in 
Part  VI.  Appendix  A contains  a derivation  of  the  fundamental 
equations  of  elastic-plastic  constitutive  models  and  is  included  for 
reference  purposes  and  future  use. 
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PART  II:  TYPICAL  MECHANICAL  BEHAVIOR  OF  SATURATED 
COHESIONLESS  SOILS 


5.  The  mechanical  response  of  a cohesionless  soil  when  subjected 
to  externally  applied  loads  is  a function  of  the  volumetric  and  devia- 
toric  stress-strain  properties  of  the  material.  These  properties,  in 
turn,  are  affected  by  such  factors  as  void  ratio,  density,  degree  of 
saturation,  interstitial  pore  fluid,  and  the  loading  history  of  the  ma- 
terial. Density,  in  particular,  has  a strong  influence  on  the  response 
of  the  material.  Loose  sand,  for  example,  contracts  and  exhibits  a 
ductile-type  stress-strain  behavior  when  subjected  to  a deviatoric  state 
of  stress,  whereas  dense  sand  dilates  and  exhibits  a brittle-type  stress- 
strain  behavior  when  subjected  to  a similar  stress  condition.  The  bound- 
ary between  these  loose  and  dense  states  is  characterized  by  that  density 
at  which  shearing  deformation  occurs  without  volume  change.  This  den- 
sity is  usually  referred  to  as  the  "critical  density,"  and,  as  shown  in 

2 

Figure  1,*  its  magnitude  varies  with  mean  normal  stress.  The 


EFFECTIVE  PRESSURE,  P' 


Figure  1.  Volumetric  response  of  sand  in  dense  and  loose  states  sub- 
jected  to  a deviatoric  state  of  stress  (constant  P') 

* For  convenience,  symbols  are  listed  and  defined  in  the  Notation 
(Appendix  B). 
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stress-strain  behavior  of  sand  is,  therefore,  highly  dependent  on  the 
confining  pressure.  For  fully  saturated  sand,  the  stress-strain  be- 
havior also  varies  greatly  depending  on  whether  the  sand  is  loaded  in  a 
drained  or  undrained  condition.  Since  the  pore  fluid  (water)  is  rela- 
tively incompressible,  pressure  builds  up  in  the  pore  fluid  during  un- 
drained loading.  Consequently,  the  effective  stress  carried  by  the  par- 
ticles comprising  the  soil  skeleton  equals  the  total  stress  minus  this 
pore  pressure.  This  effective  confining  pressure,  then,  influenciic  the 
undrained  stress-strain  behavior  of  saturated  sand. 

Behavior  Under  Monotonic-Type  Loading 

6.  Figure  2 qualitatively  depicts  typical  stress-strain-pore 
pressure  response  curves  for  four  specimens  of  saturated  sand,  each  at 

a different  initial  density,  tested  undrained  in  a triaxial  test  device. 
As  indicated  in  Figure  2,  the  specimens  are  consolidated  to  the  same 
effective  pressure  (point  l).  The  curves  marked  "l  -*■  2"  show  behavior 
typical  of  very  loose  sand.  This  specimen  develops  its  peak  strength 
at  a relatively  small  strain  (compared  with  the  value  of  strain  at  the 
end  of  the  test)  and  then  softens  (loses  strength)  with  further  strain- 
ing. The  strength  at  the  end  of  the  test,  referred  to  as  "residual 
strength,"  is  only  a small  fraction  of  the  peak  strength.  Because  of 
the  contractive  behavior  of  loose  sand,  the  effective  pressure  in  the 
specimen  decreases  as  the  test  progresses  causing  a concomitant  in- 
crease in  the  pore  pressure.  The  curves  marked  "1  ->■  5"  show  behavior 
typical  of  dense  sand.  The  strength  of  the  material  increases  contin- 
uously with  increasing  strain,  and  the  effective  pressure  in  the  dense 
sand  specimen  increases  as  the  test  is  continued,  due  to  its  dilative 
behavior.  The  pore  pressure  reaches  its  maximum  value  at  a relatively 
small  strain,  then  decreases,  eventually  becoming  less  than  its  initial 
value  (the  crossover  point  occurring  when  the  total  and  the  effective 
stress  paths  intersect). 

7.  Within  the  extreme  limits  of  very  loose  and  very  dense  sands 
a multiplicity  of  behavior  can  be  observed,  depending  on  the  initial 
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VOID  RATIO,  e „ STRESS  DIFFERENCE 


EFFECTIVE  PRESSURE,  P'  AXIAL  STRAIN, 

Figure  2.  Typical  behavior  of  saturated  sand  during  undrained  shear  in 

a triaxial  test  device 


state  of  compaction  of  the  material.  For  example,  as  indicated  by  the 
curves  marked  "1  ->■  U"  in  Figure  2,  for  a sand  slightly  looser  than  the 
densest  (dilative)  specimen,  the  effective  pressure  first  decreases  and 
then  begins  to  increase  as  the  test  continues.  Accordingly,  the  induced 
pore  pressure  first  increases  and  then  decreases  with  increasing  strain. 

8.  Typical  stress-strain  curves  for  saturated  sand  tested  under 
drained  triaxial  test  conditions  are  depicted  qualitatively  in  Figure  3. 
The  dense  sand  specimen  develops  its  peak  strength  early  in  the  test 
and  then  softens  as  the  test  progresses.  The  loose  sand  specimen,  con- 
versely, continuously  hardens  and  generally  develops  its  peak  strength 
at  the  end  of  the  test.  Comparison  of  Figure  3 with  Figure  2 indicates 
that  the  stress  difference-axial  strain  response  for  very  dense  and  very 
loose  sands  in  the  drained  condition  contrasts  markedly  with  the  corre- 
sponding responses  in  the  undrained  condition.  That  is,  a softening 
stress-strain  response  corresponds  to  a dilative  behavior  (i.e.,  dense 
sand)  in  a drained  test,  but  a contractive  behavior  (i.e.,  loose  sand) 
in  an  undrained  test.  Conversely,  a hardening  response  corresponds  to  a 
contractive  behavior  in  a drained  test  and  a dilative  behavior  in  an 
undrained  test. 


Behavior  Under  Cyclic  Loading 

9.  Typical  behavior  of  saturated  sand  subjected  to  cyclic  loading 
in  an  undrained  condition  is  shown  qualitatively  in  Figure  1+.  The  mate- 
rial is  first  consolidated  to  point  A and  then  subjected  to  small  ampli- 
tude (subfailure)  cyclic  stresses  (the  cyclic  load  oscillating  between 
zero  and  a maximum  compression  load).  The  important  feature  of  the 
behavior  depicted  in  Figure  U is  that  the  accumulated  strain  during 
cyclic  loading  prior  to  failure  (point  A -»•  B)  is  small,  compared  with 
the  strain  at  the  end  of  the  test,  whereas  the  pore  pressure  increases 
rapidly  during  this  period.  During  the  last  cycle  when  the  effective 
stress  path  reaches  the  failure  envelope  (point  B),  the  material  fails 
suddenly  and  loses  strength  as  the  test  continues  (point  B ->  C).  The 
response  of  the  materiail  from  point  B to  point  C is  similar  to  the 
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VOID  RATIO  e STRESS  DIFFERENCE 


AXIAL  STRAIN,  fj. 


Figure  3.  Typical  behavior  of  saturated  sand  during  drained  shear 

in  a triaxial  test  device 


EFFECTIVE  PRESSURE,  P ' AXIAL  STRAIN, 


Figure  h.  Typical  behavior  of  saturated  sand  under  undrained  conditions 
subjected  to  cyclic  loading  in  a triaxial  test  device 


response  of  loose  sand  under  monotonic  loading  (curves  1 -*■  2,  Figure  2); 
i.e.,  the  pore  pressure  increases  steadily  and  reaches  its  maximum  value 
at  point  C.  At  point  C the  shearing  strength  of  the  material  is  only 
a small  fraction  of  the  strength  at  point  B. 

10.  Figure  5 portrays  typical  response  of  saturated  sand  in  an 
undrained  condition  subjected  to  cyclic  stresses  with  stress  reversals 
(the  cyclic  load  oscillating  between  a maximum  extension  and  a maximxim 
compression  load).  The  most  important  feature  of  the  behavior  shown  in 
Figure  5 is  also  similar  to  that  indicated  in  Figure  U.  The  accumulated 
strain  prior  to  failure  (point  A -►  B)  is  very  small,  but  the  pore  pres- 
sure increases  rapidly  during  this  period.  At  point  B the  material 
fails  suddenly  and  loses  strength  with  fxirther  straining  while  the  pore 
pressure  increases  steadily  (point  B -*■  C). 

11.  The  behavior  of  loose  saturated  sand  in  a drained  condition 
subjected  to  cyclic  shear  stresses  with  stress  reversals  is  shown  quali- 
tatively in  Figure  6.  Tiie  rate  of  volume  change  decreases  as  the  test 
progresses;  i.e.,  most  of  the  volumetric  strain  takes  place  in  the  I’irst 
few  cycles  of  loading.  This  behavior  reflects  the  rate  of  pore  pressure 
generation  in  the  corresponding  undrained  test.  Furthermore,  the  volu- 
metric strain  in  the  first  few  cycles  is  plastic  (irrecoverable)  in 
nature  and  becomes  mostly  elastic  in  the  subsequent  cycles.  Also,  as 
indicated  in  Figure  6,  the  hysteresis  loop  generated  during  shear  is 
greatest  in  the  first  cycle  of  loading  and  becomes  smaller  as  the  test 
continues.  Such  behavior  indicates  that  during  a cyclic  drained  test, 
the  response  of  the  sand  approaches  that  of  an  elastic  material  with 
the  progression  of  the  test. 
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Figure  5.  Typical  behavior  of  saturated  sand  under  undrained  conditions 
subjected  to  cyclic  loading  with  stress  reversals  in  a triaxial  test 

device 


PART  III;  FORMULATION  OF  THE  CONSTITUTIVE  MODEL 


12.  The  mechanical  behavior  of  saturated  cohesionless  soils  was 
described  qualitatively  in  Part  II.  This  part  presents  the  development 
of  an  elastic-plastic  constitutive  relationship  that  describes  these 
various  response  features. 

13.  The  basic  theory  of  elastic-plastic  constitutive  models  for 
single-phase  (solid)  materials  is  presented  in  detail  in  Appendix  A. 
Equations  All  through  Al6  of  Appendix  A define  the  elastic  behavior  of 
the  models.  The  plastic  behavior  is  described  by  Equations  A17  through 
A28.  Equations  A29  and  A30  express  the  complete  elastic-plastic  de- 
scription. In  order  to  apply  these  equations  to  a two-phase  continuum 
(consisting  of  a solid  skeleton  and  a pore  fluid),  the  normal  stress 
components  must  be  divided  into  two  parts:  the  stress  carried  by  the 
solid  structure,  referred  to  as  the  "effective  stress,"  and  the  stress 
carried  by  the  pore  fluid,  referred  to  as  the  "pore  water  pressure." 
Mathematically,  total  stress  can  be  expressed  as 


a = o'  + u 


(la) 


or,  in  tensorieil  form,  as 


0. , 

ij 


ij 


+ u6 


iJ 


(lb) 


where 

а. .  = total  stress  tensor 
ij 

o! , = effective  stress  tensor 
ij 

u = pore  pressure 

б .  , = Kronecker  delta  = 
ij 

For  example,  in  the  case  of  a triaxial  test  (cylindrical  coordinate 
system  z , r , and  0),  Equation  lb  takes  the  following  form: 


1 i = J 
P i # j 


(2) 


lU 


T 


i 


where  and  = Og  are,  respectively,  the  axial  and  radial  total 

stress  components.  The  following  paragraphs  contain  the  mathematical 
forms  of  the  various  response  functions  developed  for  the  proposed  two- 
phase  constitutive  relation. 

Elastic  Behavior 

lU.  The  behavior  of  the  model  in  the  elastic  (recoverable)  range 
is  described  by  the  elastic  bulk  and  shear  moduli.  The  elastic  bulk 
modulus  is  assumed  to  be  a function  of  effective  mean  normal  stress,* 

P'  , or  the  first  invariant  of  the  effective  stress  tensor  (J^  = 

3P' ) (Figure  7).  The  elastic  shear  modulus,  on  the  other  hand,  is 


Figure  T.  Elastic  bulk  modulus  versus  first  invariant  of  the  effective 

stress  tensor 


assiimed  to  be  a function  of  the  second  invariant  of  the  stress  deviation 

p 

tensor,  , and  the  plastic  volumetric  strain,  (Figure  8).** 


* The  bulk  modulus  could  also  be  a function  of  the  plastic  volumetric 
strain. 

**  The  functional  forms  of  the  bulk  and  shear  moduli  (Equations  3 and 
U)  could  include  more  terms  and  provide  more  flexibility  in  fitting 
the  behavior  of  £iny  specific  material  (cf..  Reference  3). 
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Figure  8.  Elastic  shear  modulus  versus  second  invariant  of  the  stress 
deviation  tensor  showing  effects  of  plastic  volumetric  strain 


where 

K.  = initial  elastic  hulk  modulus 
1 

and  = material  constants 

G.  = initial  elastic  shear  modulus 
1 

G^,  G^,  and  G^  = material  constants 

The  constants  ^2  detemined  from  the  slope  of 

the  unloading  curve  from  an  isotropic  consolidation  test  (Figure  9)* 
The  constants  G^^  , G^  , G^  , G^  , and  Gj^  can  he  determined  from 
the  slopes  of  the  unloading  stress  difference-strain  difference  curves 
from  drained  triaxial  tests  conducted  at  different  confining  pressures 
(Figure  10 ) . 
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Figxire  9.  Proposed  relationship  for  isotropic  consolidation  test 


Plastic  Behavior 

15.  For  the  plastic  behavior,  the  loading  function  (Equa- 
tion A9)  is  assiuned  to  be  isotropic  and  to  consist  of  two  parts  (Fig- 
ure 11):  an  ui.timate  failure  envelope,  which  limits  the  maximum  shear 
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Figure  11.  Proposed  yield  surface 


strength  of  the  material,  and  a strain-hardening  yield  surface.'  The 
failure  envelope  portion  of  the  loading  function  is  assumed  to  be  of 
the  Prager-Drucker  type  and  is  denoted  by 


'(j{  j =^2  ~ “'^1 


k = 0 


and  the  strain-hardening  yield  surface  is  assumed  to  be  elliptical  and 
of  the  following  form: 


(6) 


■’i 


U<) 


2 2— 

+ R Jg  - 


X(k)  - UkT  = 0 


where  k and  a are  material  constants  representing  the  cohesive  and 
frictional  strengths  of  the  material,  respectively.  For  cohesionless 
soil  k is  negligible.  The  parameter  R is  the  ratio  of  the  major  to 
minor  axes  of  the  elliptic  yield  surface  (Figure  11 ).  Note  in  Figure  11 
that  L{k)  and  X(k)  define  the  intersections  of  each  hardening  sur- 
face with  the  failure  envelope  f , and  the  axis, 

respectively.  The  hardening  parameter  < is  generally  a function  of 
the  history  of  plastic  volumetric  strain.  For  the  present  model,  k 
was  chosen  as 


ic  = 


(7) 


.th 


he  niunber  ( 

The  following  relationships  were  chosen  between  k , L(<)  and  X(ic)  : 


where  lde,'_,_ )..  is  the  i""'  increment  of  plastic  volinnetric  strain,  and 
n is  the  niunber  of  increments  from  the  beginning  of  the  loading  history. 


L(k) 


if  ri(<)  > 0 
if  n{K)  $ 0 


(8) 


X(<)  = - i ln(l  - (9) 

When  it  is  assumed  that  the  failure  envelope  (Equation  5)  intersects 
each  ellipse  (Equation  6)  at  the  crown,  it  follows  that 


n(K) 


X(k)  - Rk 
1 + oR 


1 + oR 


(10) 
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where  D and  W are  material  constants.  W is  the  maximum  nonrecover- 
able  volumetric  compaction  that  the  material  can  withstand  under  iso- 
tropic consolidation.  Because  of  its  impact  on  the  behavior  of  the 
model  the  parameter  R is  discussed  in  the  following  paragraph. 


Parameter  R and  Its  Effect  on  Behavior  of  Model 


l6.  The  value  of  R depends  on  the  state  of  compaction  of  the 

material.  For  contractive  material  (i.e.,  loose,  Figure  2,  curve  marked 

"l  -►  2"),  the  value  of  R is  found  to  be  greater  than  1/a  . For 
dilative  material  (i.e.,  dense.  Figure  2 curves  marked  "1  -»■  U"  or 
"l  -*  5"),  the  value  of  R is  found  to  be  less  than  1/a  . R = 1/a 
corresponds  to  the  curve  marked  "l  ->  3"  in  Figure  2.  In  this  case  the 
strain-hardening  yield  surface  is  not  allowed  to  translate  relative  to 
the  origin  of  the  , axis,  and  the  behavior  of  the  model  becomes 

similar  to  the  critical  state  model  proposed  by  Schofield  and  Wroth. ^ 

These  variations  in  the  parameter  R can  be  accounted  for  by  the 
following  relation; 


R = 


R. 

1 


1 + R, 


1 + R, 


(11) 


where  ^2  material  constants  that  can  be  determined, 

for  example,  from  test  results  (through  a trial  and  error  process) 

so  that  a good  fit  to  both  the  stress  path  and  the  stress-strain  curves 
is  obtained.  Equation  11  is  used  to  compute  the  value  of  R whenever 
the  yield  surface  (Figure  ll)  is  expanding.  During  unloading  in  shear, 
however  (i.e.  path  2 ^ 3 in  Figure  12),  it  is  postulated  that  the  yield 
surface  ab  will  contract  as  shown  by  the  dashed  curve  ab'  in  Fig- 
ure 12.  It  should  be  noted  that  the  contraction  of  the  yield  surface 
is  not  dependent  on  the  plastic  strain  because,  during  unloading,  plas- 
tic strain  is  constant.  The  final  position  of  the  yield  surface  is 
determined  by  the  magnitude  of  at  point  3-  The  value  of  R at 

point  3 is  greater  than  that  given  by  Equation  11  by  a factor  R^  . 
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Figure  12.  Postulated  model  behavior  during  unloading  and  reloading 

in  shear 

Therefore  during  reloading  in  shear  (path  3 ^ in  Figure  12)  the 

parameter  R takes  the  form 

Where  the  numerical  value  of  R^  is  determined  from  the  continuity 
condition  at  the  unload /reload  interface  at  point  3.  That  is,  the  value 
of  R at  point  3 is  first  determined  from  Equations  6,  8,  and  10,  where 
X(ic)  corresponding  to  point  a is  known  and  L((c)  corresponds  to  point 
b'  . This  value  of  R , together  with  the  values  of  ^2 

can  then  be  substituted  into  Equation  12  from  which  can  be 

determined . 

IT.  The  contraction  of  the  yield  surface  as  discussed  in  the 
abo’'e  paragraph  allows  considerable  flexibility  in  modeling  the  behavior 
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oi'  strain-softening  materials  and  the  response  of  saturated  sand  under 
subfailure  cyclic  loading  conditions.  Part  IV  covers  the  behavior  of 
the  model  when  subjected  to  monotonic  and  cyclic  loadings  under  triaxial 
test . 


PART  IV:  BEHAVIOR  OF  THE  MODEL  UNDER  TRIAXIAL  TEST  CONDITIONS 


l8.  Current  methods  of  routine  stress  analysis  in  geotechnical 
engineering  require  material  property  information,  much  of  which  is 
usually  obtained  with  the  triaxial  test  apparatus.  It  is  of  interest, 
therefore,  to  investigate  the  behavior  of  the  proposed  model  in  a tri- 
axial test  configuration  for  both  drained  and  undrained  conditions, 
including  cyclic  loading.  With  the  adoption  of  the  z-axis  of  a cylin- 
drical coordinate  system  (z,  r,  and  0)  as  the  axis  of  symmetry  of  the 

assumed  isotropic  specimen,  the  effective  stress  tensor  a I and  the 

^ J 

strain  tensor  e.  associated  with  the  triaxial  test  configuration  are 
^ J 
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ij 
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(13a) 


(13b) 


The  variables  (first  invariant  of  the  effective  stress  tensor), 

Jr  (second  invariant  of  the  stress  deviation  tensor),  and  e,,  (the 
volumetric  strain  AV/V)  associated  with  the  above  stress  and  strain 
tensors  take  the  following  forms: 
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Triaxial  tests  on  saturated  sands  generally  are  performed  in  two  phases ; 
(a)  an  isotropic  consolidation  phase  and  (b)  a shear  phase.  The  shear 
phase  may  be  performed  either  drained  or  undrainea.  These  phases  are 
discussed  below. 


Isotropic  Consolidation  Phase 


19.  During  the  isotropic  consolidation  phase  of  the  triaxial  test. 
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(15a) 
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The  relation  between  the  elastic  volumetric  strain  increment  and  the 
increment  of  the  first  invariant  of  the  effective  stress  tensor  is 
(see  Equation  A12) 


dd-  = 3K  dc^j^ 


(16) 


where  the  elastic  bulk  modulus  K is  given  by  Equation  3.  Substitution 
of  Equation  3 into  Equation  I6  and  integration  of  the  resulting  expres- 
sion provides  the  following  relation  between  the  elastic  volumetric 
E 

strain  e, , and  J'  ; 

kk  1 
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SK^K. 
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(IT) 


IS 


The  relation  between  the  plastic  volumetric  strain  and 

KX 

given  by  Equation  9»  where  < (Equation  7)  for  this  phase  of  the  test 
P 


is  and  X(x)  is  , thus 
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(18) 


As  a result  of  Equations  IT  and  l8,  the  total  volumetric  strain  takes 
the  following  form; 


'kk  3K^K. 

2 1 


In 


exp(K2J|)  - 


+ W 


1 - exp 


(19) 


Equations  17  through  19  provide  a complete  specification  for  the  be- 
havior of  the  material  during  the  isotropic  consolidation  iJhase  of  the 
triaxial  test. 

20.  The  qualitative  behavior  of  the  model  during  isotropic  con- 
solidation is  shown  in  Figure  13.  The  slope  of  the  J ' versus  e,  , 

X K.K 

curve  during  initial  loading  can  be  obtained  from  Equation  19: 


Figvire  13.  Behavior  of  the  model  under  cyclic 
isotropic  consolidation 
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1 - K, 
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1 - exp 
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= 3K 


(20) 


where  K is  the  apparent  bulk  modulus  of  the  material.  The  combination 
of  Equations  3,  l8,  and  20  results  in 


K = 


^ 

i * 3D  (w  - 


4) 


(21) 


The  second  term  in  the  denominator  of  Equation  21  represents  a softening 
term;  i.e.,  it  accounts  for  the  fact  that  K is  initially  softer  than 
the  elastic  bulk  modulus  K . At  higher  pressures,  however,  this  soft- 

p 

ening  term  becomes  zero  (i.e.,  e = W)  and  the  apparent  bulk  modulus 
K equals  the  elastic  bulk  modulus  K.*  If  a sample  is  consolidated  from 
point  1 to  point  2 (Figure  13),  unloaded  from  point  2 to  point  3,  and 
then  reloaded  from  point  3 to  point  2,  the  model  dictates  that  the 
unloading-reloading  behavior  is  purely  elastic. 


Shear  Phase 


Undrained  condition 

21.  During  the  shear  phase  of  a conventional  undrained  triaxial 
test,  the  cell,  or  confining,  pressure  is  maintained  constant:  i.e., 

a = constant  = P'  (22a) 

r c 


hence. 


do  = 0 
r 


(22b) 


* For  a dense  material,  the  softening  terra  is  usually  small,  and  the 
modulus  K closely  approximates  the  elastic  bulk  modulus. 
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where  is  the  confining  pressure  at  the  end  of  the  isotropic  con- 

solidation phase.  The  stipulation  that  the  shear  phase  is  performed 
undrained  means  that  the  volumetric  strain  is  constant  and  equal  to 
the  volumetric  strain  achieved  at  the  end  of  the  consolidation  phase.* 
Accordingly,  for  an  undrained  shear  phase, 

- <*4 " '*4  ” “ 


which  leads  to 


de  = 


r 


(23b) 


The  shear  phase  up  to  failure  (Equation  5)  can  be  accomplished  by  either 
monotonic  or  cyclic  loading.  At  failure,  the  shear  strength  of  the 
material  can  either  (a)  increase  (i.e.,  undergo  strain  hardening), 

(b)  decrease  (i.e.,  undergo  strain  softening),  or  (c)  remain  constant 
(e.g.,  if  the  material  is  at  its  critical  density.  Figure  2,  point  3). 
The  behavior  of  the  model  under  these  three  conditions  is  examined 
separately  below. 

22.  Prefailure  response  under  monotonlc  loading.  The  total  volu- 
metric strain  at  the  end  of  the  consolidation  phase  of  the  test  can  be 
obtained  from  Equations  9 and  IT: 


3K2K^ 


In 


'exp^KgJj^^  - 


1 - K, 


+ w|l  - exp[-DX(K)]|  (2Ua) 


or 
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* Undrained  shear  actually  means  no  change  in  water  content;  the  con- 
dition of  constant  volumetric  strain  is  only  an  approximation. 
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The  equation  of  the  hardening  siirface  (Equation  6)  for  the  triaxial  con- 
figuration takes  the  following  form: 

2 

[j|  - L(x)]^  + |-  (o;  - - [x(k)  - L(k)]^  = 0 (25) 

An  expression  for  the  effective  stress  path  in  the  shear  phase  (i.e., 

0^  - versus  J|)  can  he  obtained  by  the  combination  of  Equations  10, 
24b,  and  25: 


+ ^(l  . cRf  0 (26) 


Figure  l4  shows  a typical  effective  stress  path  produced  by  Equation  26. 
The  sample  is  isotropically  consolidated  from  point  1 to  point  2.  The 
path  indicated  by  point  2 3 represents  the  prefailure  effective  stress 

path  under  monotonic  loading.  The  total  stress  path  is  indicated  by 
point  2^4.  Since  the  effective  and  the  total  stress  paths  are  known, 
the  excess  pore  pressure  u can  easily  be  determined  as  indicated  in 
Figure  l4.  Because  the  pore  pressure  and  the  effective  stress  tensor 
have  been  determined,  the  complete  strain  tensor  can  be  obtained  from 
Equation  A29  by  use  of  Equations  3,  4,  9»  10,  and  25  (with  R calcu- 
lated from  Equation  ll).  A computer  program  called  TDRIVER  was 
developed  to  solve  numerically  this  system  of  equations,  as  well  as 
Equation  A29,  and  to  generate  various  plots  of  stress-strain  and  pore 
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Figure  14.  Typical  effective  stress  path  produced  by  Equation  26 

pressure  response  for  undrained  triaxial  conditions.  This  program  and 
its  flow  charts  are  available  upon  request. 

23.  Postfailure  response  of  dilative  materials.  Figure  15  repre- 
sents a typical  triaxial  test  stress  path  for  dilative  material.  The 
plastic  volumetric  strain  increment  along  the  path  3 ->■  *+  can  be  obtained 
from  Equations  5,  A20,  and  A28  as 


9Ka^  + G 


Since  the  total  volumetric  strain  increment  during  the  shear  phase  is 
assumed  to  be  zero,  the  elastic  volumetric  strain  increment  can  be  set 
equal  to  the  negative  of  the  plastic  voliunetric  strain  increment;  that 


df  , -yj 

9Ka^  + G 


(hkkl)i-2  (hkk  1)2-3  ^ ( 


The  value  of  R can  be  determined  from 


R(<)  = 


Figure  15 . Typical  stress  path  for  dilative  material  under  undrained 
shear  (triaxial  test)  conditions 

The  hardening  parameter  k along  the  path  3 ^ (see  Figure  15)  is 


The  intersection  of  the  strain-hardening  yield  surface.  Equation  6,  with 
the  Jj  axis  can  be  obtained  from  Equation  9.  The  intersection  of  the 
yield  surface  with  the  failure  envelope  is  given  by 


L(ic)  = min 


X(k)  - LCk 


Note  that  when  R = 0 (i.e.,  X(k:)  = L(k)),  the  material  reaches  its 
peak  shear  strength  (point  h.  Figure  15)- 

2h.  In  summary,  for  a dilative  material,  R is  a function  of  k 
through  Equation  11,  and  k is  equal  to  the  actual  plastic  volumetric 
strain  along  the  paths  1 ->•  2 and  2 -»•  3 in  Figure  15.  Along  the  path 
3 however,  k is  given  by  Equation  29,  and  R is  a function  of  k 

through  Equation  31.  Equation  A29  will  give  the  total  strain  when  the 
computer  program  TDRIVER  is  used. 

25.  Postfailure  response  of  contractive  material.  Figure  2 
(curve  1 ->■  2)  shows  the  typical  behavior  of  a contractive  material 
under  the  undrained  triaxial  test  condition.  I'igure  I6  represents 

the  model  simulation  of  this  behavior.  The  treatment  of  path  1 •>  2 -*•  3 
is  similar  to  the  prefailure  behavior  discussed  in  the  previous  para- 
graph. Beyond  point  3 the  material  strain  softens.  This  may  be  viewed 
as  a release  of  elastic  strain  with  a concomitant  buildup  of  plastic 
strain.  The  path  3 ^ b,  therefore,  can  be  simulated  by  the  model  in  two 
steps:  elastic  unloading  from  point  3 to  point  a followed  by  elastic- 
plastic  loading  from  point  a to  point  b. 

26.  Since  path  3 -*■  a is  purely  elastic. 

Therefore,  k and  X(k:)  remain  constants  and  the  yield  surface  (Equa- 
tion 6)  is  inactive  along  this  path.  The  value  of  L , however,  is 
variable  and  can  be  computed  from  Equations  6 and  10.  The  long-dashed 
line  shown  in  the  left  sketch  of  Figure  I6  represents  the  position  of 
the  yield  surface  at  point  a.  The  value  of  R corresponding  to  this 
yield  surface  is  computed  from  Equation  12.  Along  the  path  a ->■  b,  the 
total  volumetric  strain,  which  is  constant,  can  be  obtained  from  Equa- 
tion 2l+a.  The  effective  stress  path  a -+  b can  be  determined  from  Equa- 
tion 26,  and  the  computer  program  TDRIVER  may  be  used  to  compute 
total  strain  along  the  path  3 -►  a b . 

27.  The  above  process  is  repeated  along  the  path  b ->  U as  shown 
by  the  short-dashed  line  in  Figure  16.  Note  that  the  loading  increment 
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between  3 ^ can  be  divided  into  n increments  so  that  the  dotted 

portion  of  the  stress-strain  curve  shown  in  Figure  l6  can  be  closely 
adhered  to. 


28.  Prefailure  response  under  undrained  cylic  loading.  Figures  U 
and  5 show  the  typical  behavior  of  saturated  soil  under  undrained  cylic 
shear  loading.  This  behavior  can  be  simulated  by  the  model  similarly 

to  the  manner  described  for  the  treatment  of  strain-softening  materials. 
The  path  3 a (Figiure  17)  is  purely  elastic,  in  which  ic  and  X(k) 
are  constant  while  L is  variable.  From  a to  b the  path  is  elastic- 
plastic  for  which  the  total  strain  tensor  can  be  obtained  from  Equa- 
tion A29  with  the  use  of  the  computer  program  TDRIVER. 

Drained  condition 

29.  For  the  case  of  a drained  shear  test  the  effective  stresses 
are  identical  to  the  total  stresses  and  the  stress  path  is  known  ^ priori 
Hence,  the  model  treats  this  case  identically  to  the  undrained  case 
except  that  the  restriction  of  no  volume  change  d\iring  the  shear  phase 

is  dropped.  The  material  response  is  determined  directly  from  Equation 
A29  using  the  computer  program  TDRIVER. 


PART  V:  COtJlPARISON  OF  LABORATORY  TEST  DATA  WITH  MODEL  PREDICTION 


30.  In  order  to  demonstrate  the  versatility  of  the  proposed 
constitutive  model  for  simulation  of  the  wide  variety  stress-strain-pore 
pressure  responses  observed  in  tests  of  saturated  cohesionless  soils  (as 
seen  in  Figures  2 and  I*),  four  sets  of  experimental  data  were  considered. 

TI»e  first  set  of  data  (data  set  l)  typifies  the  behavior  of  contractive 
(loose)  materials.  As  indicated  in  Figure  2 (curves  1 -*•  3),  this  type 
of  material  exhibits  an  increase  in  pore  pressure  during  undrained  shear. 

The  second  set  of  data  (data  set  2)  typifies  the  behavior  of  dilative 
(dense)  materials.  During  undrained  shear  these  materials  exhibit  an 
increase  in  effective  pressure  and  a consequent  decrease  in  pore  pres- 
sure (curves  1 ->  5,  Figure  2).  The  third  set  of  data  (data  set  3) 
typifies  the  behavior  of  very  loose  (work-softening)  materials  subjected 
to  undrained  shear  (curves  1 -*•  2,  Figure  2).  The  fourth  set  of  data 
(data  set  I4 ) demonstrates  the  behavior  of  loose  saturated  sand  under 
undrained  cyclic  shear  loading.  As  shown  in  Figure  U,  undrained  cyclic 
shear  causes  a rapid  decrease  in  effective  pressure  and  a consequent 
increase  in  pore  pressure.  The  model  can  also  simulate  the  response  of 
dense  saturated  sand  under  undrained  cyclic  shear  loading.  However, 
because  of  lack  of  experimental  data  in  proper  format,  such  simulation 
was  not  made. 

Material  Constants 

31-  The  numerical  values  of  the  material  constants  used  for  model 
simulation  of  the  four  sets  of  data  are  given  in  Table  1.  As  indicated 
in  this  table,  there  are  15  material  constants  associated  with  the 
constitutive  model.  Three  material  constants  (K^  , , and  K^)  are 

associated  with  the  elastic  bulk  modulus  (Equation  3),  and  five  mate- 
rial constants  (G^  , , G^  , G^  , and  Gj^)  are  associated  with  the 

elastic  shear  modulus  (Equation  k) . The  material  constants  k and  a 
relate  to  the  ultimate  failure  envelope  of  the  material  (Equation  5)» 
and  the  material  constants  W and  D are  associated  with  the  hardening 
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function  (Equation  9).  The  material  constants  R.  , R.  , and  R are  ' 

1 i <;  ^ 

associated  with  the  parameter  R (Equation  11),  which  defines  the  ratio 
of  the  major  to  the  minor  axes  of  the  elliptic  yield  surface.  Direct 

determination  of  these  material  constants  requires  an  elaborate  set  of  ^ 

experimental  data.  Because  of  the  lack  of  detailed  experimental  data 
some  of  the  material  constants  in  Table  1 were  assumed  to  be  zero.  In 

addition,  the  numerical  values  of  some  of  the  other  constants  had  to  be  ; 

determined  by  trial  and  error,  by  the  use  of  the  computer  program 
TDRIVER  in  order  to  obtain  a good  fit  to  the  available  experimental  data, 
in  lieu  of  the  values  being  determined  directly  from  test  results.  Ac- 

cordingly,  the  entire  capability  of  the  constitutive  model  was  not  1 

utilized  in  simulating  these  experimental  data.  The  following  I 

sections  compare  the  model  predictions  with  these  experimental  data.  ^ 

( 

Comparison  with  Contractive  Material  (Data  Set  l)  , 

32.  The  experimental  data  for  the  contractive  material  (Reid- 
Bedford  Model  sand)  were  obtained  from  Reference  ^4.  The  data  consist 
of  results  from  one  load/unload  isotropic  consolidation  test  and  two 
consolidated  undrained  triaxial  tests.  The  tests  were  conducted  on 
saturated  samples  of  Reid-Bedford  Model  sand  having  an  initial  relative 
density  of  approximately  20  percent.  Figures  l8  through  21  compare  the 
actual  test  results  with  predicted  behavior  by  the  model  for  this  set  of 
data.  The  result  of  the  consolidation  test  is  shown  in  Figure  l8.  Fig- 
ure 19  shows  effective  stress  paths  in  principal  stress  difference- 
effective  mean  normal  stress  space.  The  shear  stress-strain  relations, 
in  terms  of  principal  stress  difference  versus  axial  strain,  are  shown 
in  Figure  20.  Figure  21  shows  excess  pore  pressures  versus  axial  strain. 

Figures  l8  through  21  indicate  that  the  constitutive  model  qualitatively 
simulates  the  response  of  contractive  materials. 

Comparison  with  Dilative  Material  (Data  Set  2) 

33.  Experimental  data  for  a dilative  material  were  obtained  from 
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Figure  l8.  Isotropic  consolidation  resKnse;  experi 
mental  versus  model;  data  set  1 


Figure  19.  Comparison  of  effective  stress  paths 
for  consolidated  undrained  triaxial  tests;  experi- 
mental versus  model;  data  set  1 


EXCESS  PORE  PRESSURE 


Reference  ,1.  The  data  consist  of  the  same  type  of  test  results  as  re- 
ported for  the  contractive  material.  The  tests  were  conducted  on  satu- 
rated samples  of  Reid-Bedford  Model  sand  having  an  initial  relative 
density  of  approximately  T6  percent.  Figures  22  through  25  compare  the 
actual  test  results  with  behavior  predicted  by  the  model  for  the  dila- 
tive material.  The  result  of  the  consolidation  test  is  shown  in  Fig- 
ure 22.  The  effective  stress  paths  are  presented  in  Figure  23,  and 
Figure  2h  shows  principal  stress  difference  versus  axial  strain  rela- 
tions. Figure  25  shows  excess  pore  pressures  versus  axial  strain.  In 
generating  these  model  simulations,  the  yield  function  F (the  elliptic 
cap)  was  used  as  the  loading  function  from  the  beginning  of  the  test 
until  the  effective  stress  path  reached  the  ultimate  failure  envelope 
f ; from  then  on,  the  ultimnte  failure  envelope  was  used  as  the  loading 
function.  During  the  first  part  of  the  model  simulation  (i.e.,  when  F 
was  used  as  the  loading  function)  pore  pressure  increased  and  reached 
its  maximum  value  when  the  effective  stress  path  reached  the  ultimate 
failure  envelope.  As  the  test  continued,  excess  pore  pressure  then 
decreased  and  eventually  became  negative.  The  experimental/model  com- 
parisons shown  in  Figures  22  through  25  indicate  that  the  proposed 
constitutive  relation  can  qualitatively  simulate  the  stress-strein- 
pore  pressure  response  of  dilative  materials. 

Comparison  with  Strain-Softening  Behavior  (Data  Set  3) 

3^.  Experimental  data  for  a strain-softening  material  were  ob- 
tained from  Reference  5.  The  data  consist  of  results  obtained  from  one 
consolidated  undrained  triaxial  test.  The  test  was  conducted  on  a satu- 
rated sample  of  Banding  sand  having  an  initial  relative  density  of  ap- 
proximately 13  percent.  Comparisons  of  the  test  results  with  behavior 
predicted  by  the  model  are  shown  in  Figures  26  through  28.  Figure  26 
shows  the  effective  stress  paths,  and  Figure  27  presents  principal 
stress  difference  versus  axial  strain  relations.  Excess  pore  pressure 
versus  axial  strain  is  shown  in  Figure  28.  These  figures  illustrate 
that  the  proposed  constitutive  relation  quite  capably  simulates  the 
stress-strain-pore  pressure  response  of  work-softening  materials. 
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MODEL  BEHAVIOR 


Figure  27.  Principal  stress  difference  versus  axial  strain 
response  from  undrained  triaxial  test;  experimental  versus 

model;  data  set  3 


INITIAL  EFFECTIVE  CONFINING 
PRESSURE  = 4.  0 kg/cm2 


Figrire  28.  Excess  pore  water  pressure  versus 
axial  strain  response  from  undrained  triaxial 
test;  experimental  versus  model;  data  set  3 

Comparison  with  Cyclic  Test  Results  (Data  Set  4) 


35.  Cyclic  test  data  were  obtained  from  Reference  5 and  consist 
of  results  from  one  consolidated  undrained  triaxial  test.  The  test  was 
conducted  on  a saturated  sample  of  Banding  sand  having  an  initial  rela- 
tive density  of  approximately  21  percent.  The  cyclic  load  oscillated 
, between  zero  and  a maximum  compression  load.  Figures  29  through  31 

show  compeirisons  of  the  test  results  with  behavior  predicted  by  the 
model  for  this  set  of  data.  Figure  29  shows  the  effective  stress  paths. 
The  principal  stress  difference  versus  axial  strain  responses  are  shown 
in  Figure  30,  and  excess  pore  pressure  versus  axial  strain  is  compared 
( in  Figiire  31.  The  most  significant  experimental  feature  of  cyclic  load- 


ing  is  the  progressive  increase  of  pore  pressure  observed  with  increas- 
ing stress  cycles , a phenomenon  which  contributes  to  liquefaction  of 


I paths  for  consolidated 

undrained  (cyclic  shear)  triaxial  test;  experimental 
versus  model;  data  set  1* 


EXCESS  PORE  PRESSURE  a.  kg/cm2 


saturated  sand.  Figure  31  shows  that  the  proposed  constitutive  relation 
simulates  this  phenomenon  quite  well. 


Figure  31.  Excess  pore  water  pressure  versus  axial  strain 
response  from  cycled  undrained  triaxial  test;  experimental 
versus  model;  data  set  h 
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PART  VI;  SUMMARY  AND  RECOMMENDATIONS 


36.  A three-dimensional  elastic-plastic  isotropic  constitutive 
relation  has  been  developed  for  simulating  the  stress-strain-pore 
pressure  response  of  saturated  cohesionless  materials.  The  major  new 
features  of  the  constitutive  model  are  (a)  its  treatment  of  work- 
softening behavior,  and  (b)  its  ability  to  predict  progressive  increases 
of  pore  pressure  under  low-amplitude  (subfailure)  cyclic  shear  loadings. 
Both  of  these  phenomena  are  important  factors  that  contribute  to  the 
liquefaction  of  saturated  sands.  The  constitutive  relation  also  allows 
for  dependency  of  the  elastic  shear  modulus  on  the  shearing  stress  and 
plastic  volumetric  strain.  This  particular  feature  of  the  model  is 
quite  useful  for  simulating  the  nonlinear  shearing  stress-strain 
response  of  soils. 

37.  The  formulation  of  the  work-softening  behavior  is  numerically 
stable.  Theoretically,  however,  this  behavior  could  violate  the 
Drucker^  stability  postulate,  which  is  a sufficient,  though  not  neces- 
sary, condition  for  satisfying  all  of  the  thermodynamic  and  continuity 
requirements  of  the  incremental  theory  of  plasticity. 

38.  It  is  recommended  that  Phase  III  of  this  investigation 
include  incorporation  of  this  constitutive  model  into  a suitable 
numerical  computer  code  for  use  in  the  assessing  of  the  liquefaction 
potential  of  dams  and  other  earth  structures  subjected  to  different 
transient  or  static-type  loadings. 
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APPENDIX  A:  FUNDAI4ENTAL  EQUATIONS  OF  ELASTIC-PLASTIC 

MATERIAL  MODELS  - 

Basic  Concepts  from  Continuum  Mechanics 

1.  In  engineering  practice  it  is  convenient,  and  often  reasonable, 
to  disregard  the  structural  details  of  materials  and  consider  their 
gross  behavior  only.  Engineering  materials  are  therefore  described,  or 
characterized,  mathematically  within  the  frameworks  of  the  theory  of 
continuous  mass  media.  Neglecting  thermal  effects,  the  basic  field 

equations  that  govern  the  motion  of  a continuum  are  the  continuity  ; 

equation* 


iP 

9t 


+ 


(Al) 


and  the  equations  of  motion 


F.  - pa.  = 0 
1 1 


(A2) 


where 


p = mass  density 
t = time 

v^  = components  of  velocity  vector 
= symmetrical  stress  tensor 
F^  = components  of  body  force 
a^  = components  of  acceleration  vector 
2.  Equations  Al  and  A2  constitute  four  equations  that  involve  ten 
unknown  functions  of  time  and  space;  The  mass  density  p , the  three 


velocity  components 


v^  , and  the  six  independent  stress  components 


o. 


ij 


The  body  force  components  F^  are  known  quantities  and  the 


acceleration  components  a.  are  expressible  in  terms  of  the  velocity 


* Indices  assume  values  1,  2,  or  3.  A repeated  index  is  to  be  summed 
over  its  range.  A comma  in  the  subscripts  represents  a derivative. 
Quantities  are  referred  to  rectajigular  Cartesian  coordinates  . 
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components  . Therefore,  six  additional  equations  relating  the  cen 
unknown  variables  are  required  in  order  to  determine  the  notion  or 
deformation  of  a medium  when  subjected  to  external  disturbances  such 
as  surface  forces  and/or  displacements.  In  continuum  mechanics,  such 
relations  are  stated  by  constitutive  equations  (or  material  models)  which 
relate  stresses  to  deformation  and  history  of  deformation.  The  dif- 
ference between  constitutive  equations  and  field  equations  (Equations  A1 
and  A2)  is  that  the  latter  are  applicable  to  all  materials,  whereas  the 
former  represent  the  intrinsic  response  of  a particular  material  or 
class  of  materials. 

3.  The  general  form  of  a constitutive  equation  may  be  expressed 
by  the  functional  form 
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where  the  deformation-rate  and  spin  tensors,  D and  SI  , respec- 

mn  qp 

tively,  are  related  to  the  components  of  the  velocity  vector  v^ 
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and  the  strain  tensor  e is  related  to  the  components  of  the  displace- 
rs 

ment  vector  - For  small  displacement  gradients. 
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Equations  A1  through  A3  constitute  ten  equations  which  include  ten 
unknown  variables.  These  equations  will  lead,  in  conjunction  with  the 
kinematic  relations  given  by  Equations  A.U  and  A5  and  boundary  conditions, 

to  a complete  description  for  solution  of  a boundary-value  problem. 

I4.  The  mechanical  behavior  of  a number  of  engineering  materials 
is  described  within  the  framework  of  elastic-plastic  constitutive 
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relationships.  The  development  of  the  specific  functional  form  of 
Equation  A3  for  the  elastic-plastic  models  is  given  in  the  following 
section. 


General  Description  of  Elastic-Plastic  Constitutive  Models 

5.  The  basic  premise  of  elastic-plastic  constitutive  models  is 
the  assumption  that  certain  materials  are  capable  of  undergoing  small 
plastic  (permanent)  as  well  as  elastic  (recoverable)  strains  at  each 
loading  increment.  Mathematically,  the  total  strain  increment  is 
assumed  to  be  the  sum  of  the  elastic  and  plastic  strain  increments;  i.e.. 


de . , = de . , + de . , 
ij  ij  ij 


(a6) 


where 

de^j  = components  of  the  total  strain  increment  tensor 
E 

de . . = components  of  the  elastic  strain  increment  tensor 
^ J 

p 

de^j  = components  of  the  plastic  strain  increment  tensor 

6.  Within  the  elastic  range  the  behavior  of  the  material  can  be 
described  by  an  elastic  constitutive  relation  of  the  type 


de?,  = A.  ,,  , (o  ) do, - 
ij  ijkl  \ mn/  kl 


(A7) 


where 


*ijkl 


do 


kl 


material  response  function 
components  of  stress  increment  tensor 


The  behavior  of  the  material  in  the  plastic  range  can  be  described 

within  the  framework  of  the  generalized  incremental  theory  of  plasticity. 

6* 

The  mathematical  basis  of  the  theory  was  established  by  Drucker  by 


* Raised  numerals  refer  to  similarly  numbered  items  in  the  References 
at  the  end  of  the  main  text. 
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introducing  the  concept  of  material  stability  with  the  following  impli- 
cations : 


a.  Yield  surface  (loading  function)  should  be  convex  in 
stress  space. 

b.  Yield  surface  and  plastic  potential  should  coincide  (this 
is  referred  to  as  associated  flow  rule). 

c_.  Work  "softening"  should  not  occur. 

These  three  conditions  can  be  summarized  mathematically  by  the  following 
inequality: 


da. , dE . . > 0 
ij  ij  - 


(A8) 


The  above  conditions  allow  considerable  flexibility  in  choosing  the  form 
of  the  loading  function  for  the  model,  which  serves  as  both  a yield 
surface  and  plastic  potential.  For  isotropic  materials  the  yield  sur- 
face may  be  expressed,  for  example,  as 


(j,  ,vr, , «) 


(A9) 


where 


J,  = o = first  invariant  of  stress  tensor 
1 nn 

Jg  = ^ = second  invariant  of  stress  deviation  tensor 

S.  = a.,  - (J, /3)6..  = stress  deviation  tensor 
ij  ij 


Kronecker  delta 


K = a hardening  parameter 

The  hardening  parameter  < , in  general,  can  be  taken  to  be  a function 

p 

of  the  plastic  strain  tensor  e^j  . The  yield  surface  of  Equation  A9 
may  expand  or  contract  as  the  hardening  parameter  tc  increases  or 
decreases,  respectively  (Figure  Al). 

7.  Conditions  a,  b,  and  c above,  taken  in  conjunction  with 
Equation  A9,  result  in  the  following  plastic  flow  rule  for  isotropic 
materials : 


Al* 
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Figure  Al.  Typical  yield  surface  for  an  elastic- 
plastic  model 
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ij 
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where  dX  is  a positive  scalar  factor  of  proportionality,  which  is  non- 
zero only  when  plastic  deformations  occur,  and  is  dependent  on  the  par- 
ticular form  of  the  loading  function. 

Elastic  strain  increment  tensor 

8.  For  isotropic  elastic  materials  the  strain  increment  tensor 
(Equation  A7)  takes  the  following  form: 


(All) 


The  bulk  and  shear  moduli  can  be  functions  of  the  invariants  of  the 
stress  tensor.  Accordingly,  assume  that  K = K and 

G = » *^2  ’ *^3)  ’ *^3  third  invariant  of  the  stress 


deviation  tensor.  Equation  All  can  be  written  in  terms  of  the  hydro- 
static and  deviatoric  components  of  the  strain  and  stress  increment 
tensors ; i . e. , 


where 

E 

de,  , = increment  of  elastic  volumetric  strain 
kk 

E 

deT.  = elastic  strain  deviation  increment  tensor 
ij 

In  order  not  to  generate  energy  or  hysteresis  within  the  elastic  range, 
the  elastic  behavior  of  the  model  must  be  path  independent.  The  mate- 
rial should  then  possess  a positive  definite  elastic  internal  energy 
function  f , which  is  independent  of  stress  path.  The  strain  energy 
function  can  be  written  as 


In  order  for  4*  to  be  independent  of  stress  path,  the  integrals  in 
Equation  AlU  have  to  depend  only  on  the  current  values  of  J,  and  Jp  . , 


Therefore,  the  bulk  and  shear  moduli  have  to  be  expressed  as 


K = K 


G = G 


(A15) 


Further,  K and  G must  always  be  positive.  Since  during  elastic 
deformation  the  hardening  parameter  < is  constant,  the  bulk  and  shear 
moduli  can  also  be  expressed  as 


(Al6) 


Plastic  strain  increment  tensor 

9.  The  plastic  strain  increment  tensor  is  given  by  Equation  AlO 

where  the  loading  function  is  given  by  Equation  A9.  The  hardening 

function  in  Equation  A9  could  be  taken  as  being  equal  to  plastic  volu- 
p 

metric  strain  ; thus 


K = 


(AI7) 


The  use  of  Equation  AIT  will  allow  the  yield  surface  to  expand  as  well 
as  to  contract.  Figure  Al.  The  plastic  loading  criteria  for  the 
function  are  given  as 


> 0 loading 
= 0 neutral  loading 
< 0 unloading 


(AI8) 


Because  deT,  = 0 during  unloading  or  neutral  loading,  as  well  as  for 
< 0 , Equations  All  through  A13  are  used  to  determine  the  purely  elas- 
tic strain  changes.  The  prescription  that  neutral  loading  produces  no 
plastic  strain  is  called  the  continuity  condition.  Its  satisfaction 


leads  to  concidence  of  elastic  and  plastic  constitutive  laws  during 
6 T 

neutral  loading.  ’ 

10.  Like  the  elastic  behavior,  the  plastic  stress-strain  relation 
can  be  expressed  in  terms  of  the  hydrostatic  and  deviatoric  components 
of  strain.  Application  of  the  chain  rule  of  differentiation  to  the 
right  side  of  Equation  AlO  yields 
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Multiplication  of  both  sides  of  Equation  A19  by  6. . gives 


de^  = 3 dX  1^ 
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The  deviatoric  component  of  the  plastic  strain  increment  tensor  (de..) 
can  be  written  as 


(A21) 


Substitution  of  Equation  A19  and  A20  into  Equation  A21  yields 


= -~=  S.  , (A22) 
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11.  In  order  to  use  Equations  A20  and  A22,  or  Equation  AI9,  the 
proportionality  factor  dX  must  be  determined.  This  can  be  accomplished 
in  the  following  manner.  From  Equations  A9  and  A17  the  total  derivative 
of  (5  becomes 
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In  view  of  Equations  A12,  A13,  and  A20,  Equation  A23  "becomes 
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Substitution  of  Equation  A6  into  Equation  A2h  results  in 
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With  the  substitution  of  the  values  of  de^^j^  and  de^j  from  Equa- 
tions A20  and  A22,  respectively.  Equation  A26  becomes 
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When  the  equation  is  solved  fvor  dX  , it  is  expressed  as* 
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» For  elastic-ideally  plastic  material,  the  numerator  in  Equation  A28 
Is  d({  and  * 
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Total  strain  increment  tensor 


12.  The  total  strain  increment  tensor  can  be  obtained  by  combin- 
ation of  Equations  a6,  All,  A19,  and  A28;  thus. 
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Similarly,  the  stress  increment  tensor  can  be  written  as 
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Equation  A29,  or  Equation  A30,  is  the  general  constitutive  equation  for 
an  elastic-plastic  isotropic  material.  To  use  these  equations  it  is 
only  necessary  to  specify  the  functional  forms  of  K , G , and  (5  . 
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APPENDIX  B:  NOTATION 
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Components  of  acceleration  vector 
Material  response  function 

Flastic  strain  deviation  increment  tensor 
Components  of  the  total  strain  increment  tensor 

Components  of  the  elastic  strain  increment  tensor 

Increment  of  elastic  voliometric  strain 

Components  of  the  plastic  strain  increment  tensor 

Increment  of  plastic  volumetric  strain 

Positive  scalar  factor  of  proportionality  used 
with  flow  rule 

Components  of  stress  increment  tensor 
Material  constant 
Deformation  rate  tensor 
Void  ratio 

Failure  envelope  portion  of  the  loading  function 

Loading  function 

Components  of  body  force 

Elastic  shear  modulus 

Initial  elastic  shear  modulus 

Material  constants 

First  invariant  of  effective  stress  tensor 
Second  invarieint  of  the  stress  deviation  tensor 

Third  invariant  of  the  stress  deviation  tensor 

Materieil  constant  representing  cohesive  strength 
of  material 

Elastic  bulk  modulus 
Apparent  bulk  modulus  of  material 
InitieLL  elastic  btilk  modulus 
Material  constants 

Effective  pressure  or  effective  mean  normaJ. 
stress 
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R.  ,R. ,R„ 
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Confining  pressure  at  end  of  isotropic  consoli- 
dation phase 

Ratio  of  the  major  to  the  minor  axis  of  the 
elliptic  yield  surface 

Change  in  R during  cyclic  loading 

Material  constants 

Stress  deviation  tensor 


u Excess  pore  pressure 

u^  Components  of  displacement  vector 

v^  Components  of  velocity  vector 

W Material  constant,  maximum  nonrecoverable 

volumetric  compaction  material  can  withstand 
under  isotropic  consolidation 

z,r,0  Cylindrical  coordinate  system 

a Material  constant  representing  frictional 
strength  of  material 

AV/V  Volumetric  strain 

6 . . Kronecker  delta 

ij 

e. . Strain  tensor 

ij 

Volximetric  strain 

e Radial  strain 

r 

e Axial  strain 

z 

Gg  Tangential  strain  components 

p 

e, , Plastic  volumetric  strain 

kk 

K Hardening  parameter 

p Mass  density 

a..  Total  stress  tensor 

ij 

a!.  Effective  stress  tensor 

ij 

Radial  total  stress  component 

Axial  stress  component 

a„  Tangential  total  stress  component 
0 

V Complementary  energy  density  function 

^ Spin  tensor 
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